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We present a new method for the calculation of frag- 
ment size correlations in a discrete finite system in which 
correlations explicitly due to the finite extent of the sys- 
tem are suppressed. To this end, we introduce a com- 
binatorial model, which describes the fragmentation of a 
finite system as a sequence of independent random emis- 
sions of fragments. The sequence is accepted when the 
sum of the sizes is equal to the total size. The param- 
eters of the model, which may be used to calculate all 
partition probabilities, are the intrinsic probabilities asso- 
ciated with the fragments. Any fragment size correlation 
function can be built by calculating the ratio between the 
partition probabilities in the data sample (resulting from 
an experiment or from a Monte Carlo simulation) and the 
'independent emission' model partition probabilities. This 
technique is applied to charge correlations introduced by 
Moretto and collaborators. It is shown that the percola- 
tion and the nuclear statistical multifragmentaion model 
(smm) are almost independent emission models whereas 
the nuclear spinodal decomposition model (bob) shows 
strong correlations corresponding to the break-up of the 
hot dilute nucleus into nearly equal size fragments. 

PACS numbers: 25.70.Pq, 24.60.Ky, 05.10.-a, 24.10.Lx 



I. INTRODUCTION 

The break-up of any finite composite sys- 
tem (atomic clusters, atomic nuclei, fullerenes, 
molecules. . . ) is characterized by a probability dis- 
tribution which incorporates constraints imposed by 
dynamical or static conservation laws. Thus, in the 
case of nuclear decay, the observed multifragmenta- 
tion modes provide information on properties of nu- 
clear matter at high excitation energy. From a statis- 
tical view point the simplest fragmentation model may 
be formulated by attributing an independent emission 
probability to each type of fragment (mass, charge). 
In the limit of infinite parent system size the resulting 
model (which will be referred to herein as the indepen- 
dent emission model) exhibits no correlation between 
fragments. For finite systems, we show hereafter that 
the correlations induced by the static conservation 



laws, that is in mass and/or charge)] (referred to here- 
after as trivial correlations), can be exactly calculated. 
In the independent emission model, all the physical in- 
formation is contained in the emission probabilities of 
the different types of particles. However, because of 
the static conservation laws, these intrinsic probabili- 
ties are not equal to the observed probabilities. 

Most theoretical multifragmentation models, which 
describe the process of instantaneous break-up of the 
atomic nucleus submitted to extreme temperature and 
pressure conditions, introduce other forms of corre- 
lations between particle types. When these correla- 
tions are specific to a given model, their experimen- 
tal observation constitutes a crucial test of valida- 
tion/invalidation. For example, several models, de- 
scribe the decay of hot nuclei by the development 
of density fluctuations (surface or volume instabilities 
jl] [I)). Among these models, the spinodal nuclear de- 
composition mechanism may dominate when the col- 
lision between two nuclei leads to a highly heated 
and sufficiently compressed nucleonic system. The 
decompression phase leads the system into the spin- 
odal zone (zone in which the incompressibility modu- 
lus (dP/dp)s is negative) where the density fluctua- 
tions are exponentially amplified up to produce frag- 
mentation |5J. The dynamics of the density waves in 
the system is dominated by the most un-stable mode, 
whose wave length is of the order of 10 fm. Thus the 
composite nucleus will disintegrate into almost equal 
size fragments (in the range Z « 6 — 20 and, more 
particularly, Z w 10 — 15, see Section IV D). 

Binary sequential de-excitation models (gemini Q 
or SIMON |7|,|8|]) do not exhibit, of course, any pref- 
erential decay into equal charges. Nor do instanta- 
neous multifragmentation models : the Copenhagen- 
Moscow model (code smm) || (see Ref. ju| for t he 
system Xe+Sn at 32 A MeV, and the Secti on |IV B| of 
the present article) and the Berlin model |ll] (code 

MMMC) ]T^| . 

Experimentally, the charge distributions are privi- 
leged tools for the the study of nuclear multifragmen- 
tation. However, the yields of various charges alone 



*In practice, because of the difficulties with mass mea- 
surements, studies are mainly carried out on charge par- 
titions, noted n : (ru , . . . , n ztot ) where n z is the number 
of charges z in the partition. The charge conservation law 
reads : z Uz = Ztot • 
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do not permit a sufficient discrimination of mecha- 
nisms. Model validations thus require the comparison 
of intra-event charge correlations. In this context, a 
difficulty arises from the fact that the detected frag- 
ments are not produced only at the multifragmen- 
tation stage of the reaction. Certain light particles, 
for example, are emitted during the inter-penetration 
of the nuclear spheres (pre-equilibrium phase), oth- 
ers are emitted, at the end of the process, by the hot 
multifragmentation fragments. The final (detected) 
partitions have thus, in part, lost the memory of the 
crucial moment of the reaction. Therefore, it may 
be necessary to use statistical methods in order to 
detect the charge correlations induced by the initial 
multifragmentation. One of these methods, proposed 
by L. Moretto and collaborators |U5), was shown to 
be especially efficient for detecting the presence of 
the spinodal decay mechanism (volume instabilities) 
of the nucleus. 

The fragments formed during the spinodal decom- 
position phase have comparable sizes (charges) . How- 
ever, this effect is not visible in the charge spectra gen- 
erated by a Monte-Carlo code (Brownian One Body 
Dynamics, BOB code [m) simulating this type of mech- 
anism. The reasons are numerous : coalescence and 
primary fragment deexcitation, finite size effect induc- 
ing mode superpositions . . . The same remark applies 
to the experimental charge distributions : no excess 
yield is visible in the expected charge domain. How- 
ever, the method of charge correlations reveals, for 
this model, a small "fossil" signal which corresponds 
to events in which the system breaks into similar size 
fragments and whose charges have not been modified 
(or reduced by the same quantity) before detection. 
The method introduced by Moretto and collaborators 
consists in calculating the correlation function of the 
mean charge (Z) of the M imf[] and of their standard 
deviation a. A peak appears therefore in this corre- 
lation function for a « and (Z) « 10 — 15. Experi- 
mentally, a peak has been effectively observed for the 
Xe+Sn system at 32 A MeV in central collisions with 
the INDRA multidetector |I(],[l4|]. Preferential decom- 
positions in three approximately equal size fragments 
were also observed in central Xe+Cu reactions at 45 
A MeV with the multics multi-detector |l5) . 

The goals of this article are the following : 

• We wish to make the interpretation of correla- 
tions more rigorous. Progress is necessary be- 
cause the peak related to spinodal decomposi- 
tion is often generated by a very small number 



intermediate Mass Fragments. Concretely fragments 
with charge greater or equal to a given limit (Zmi n = 3 
or 5 in this work). 



of experimental or synthesized events. It corre- 
sponds, as we will see, to the ratio of two very 
small quantities and therefore will be character- 
ized by a large error bar. The significance as- 
sociated with a peak must therefore be system- 
atically evaluated. To this end, we show that 
the error in the denominator of the correlation 
function can be greatly reduced by substituting 
a convolution product for the random selection 
process proposed in the initial method. 

• The correlation peak corresponding to the spin- 
odal decomposition (or to any other cause) is su- 
perimposed on a dominant structure due to the 
correlations induced by the total charge conser- 
vation law (trivial correlations). This structure 
often makes the interpretation of the peaks in 
terms of physically interesting correlations diffi- 
cult or ambiguous. Hence, it is important to 
correct the correlation function for finite size 
effects. For this reason, it has been proposed 
to construct the denominator of the correla- 
tion function in a different way than that intro- 
duced in Ref. Jl3|] using the minimum informa- 
tion model. It will be shown that this method 
can hide peaks corresponding to non trivial cor- 
relations. 

• We therefore introduce, in an algebraical exact 
way, the effects of charge conservation using the 
independent emission model. Thanks to this 
new method, any event sample with only trivial 
correlations will show a flat correlation function. 

• We will study more completely the independent 
emission model constrained by the charge con- 
servation. The notion of intrinsic probability of 
particles will be introduced. 

• Finally, this new method will be validated by 
its application to three nuclear decomposition 
models (SMM, percolation |l6| and bob). It will 
be shown that these models are, to first order, 
independent emission models. 

II. THE CHARGE CORRELATION FUNCTION 

A. Algebraic calculation of the denominator 

The quan- 
tity l + R{(Z),a\M) = P({Z),a\M)/P nc ({Z),a\M), 
where (Z) is the mean charge of the IMF, a their 
standard deviation and M their multiplicity, will be 
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called the charge correlatior^ The method tradition- 
ally used to calculate the denominator of a correlation 
function [|l3| consists in constructing "pseudo-events" 
using randomly selected fragments belonging to dif- 
ferent events of the sample with a given IMF multi- 
plicity. The global variable distributions relative to 
the pseudo-events do not contain intra-cvent corre- 
lations. The numerator P and the denominator P uc 
of the correlation function are calculated in the same 
way, the first one from the sample events, the sec- 
ond one from the pseudo-events. Since the denom- 
inator does not contain intra-event correlations, its 
probability density function is written with an index 
uc (un-correlated). 

The only experimental information required for the 
calculation of the denominator is the charge distri- 
bution of the sample. It is equivalent, and, from a 
computational point of view, faster, to sort charges 
with respect to the average charge distribution, rather 
than to select fragments among events. In fact, the 
random selection using the charge probability distri- 
bution is not even necessary since the denominator 
can be calculated algebraically in the form of a con- 
volution product. 

One notes P((Z)\M) the probability to obtain a 
value (Z) of the mean IMF charge for the multiplicity 
M events (z2(z) p (( z ) \ M ) = 1 > hereafter all condi- 
tional probabilities will be assumed to be normalized 
by a relation of the same type). This conditional prob- 
ability is given by the convolution : 

P UC «Z)|M) = £... J2 Pz{Zi\M) 

. ..p z (z m -i\m) p z {m(z) - (m - i)(zy \M), (1) 

where (Z)' is the mean charge of the IMF except 
the last and P Z (Z\M) the IMF charge probability 
distribution for a given multiplicity. The last fac- 
tor accounts for the IMF total charge conservation 

-,M 



ability to obtain a standard deviation a, when the 
fragments are randomly selected, is : 

P uc (a\M) = J2--- E Pz{Zx\M)... 

Zx 

xP z (Z M -i\M) P z {(Z)' + M 
x P Z ({Z)'-M 



r'2 



M-l M 



a 2 a' 2 

-1 _ M) , (3) 

M —I M ' ' w 



where a' is the standard deviation of the charges of 
the IMF except the last. If the term under the square 
root is negative, the probability is zero. Finally the 
correlation between the mean charge and the standard 
deviation reads : 



P uc ((Z),a \M) = 



Zm 



Pz{Zi\M). 



xPz(Zm-i\M) P z (Z m \M) 6 Zu ,m(z)+(m-x)(z)' 



Z M ,{Z)'±M yf- 



(4) 



where 5 a ,b, the Kronecker symbol, is equal to 1 when 
a = b and otherwise. The multinomial decomposi- 
tion leads to an equivalent (but more practical) form 
of this equation : 



E 



P uc ((Z),a \M) = 
P z (Z\M) Nz 



n 



N z l 



(5) 



Ez N z = M 

^2, z ZN Z = M{Z) 

J2z z2N z = M i(z} 2 +^ 2 ) 

where N z is the number of IMF with charge Z and 
N an IMF partition. The product runs over all possible 
IMF charges. The probabilities in the denominator re- 
spect the normalization : J2aJ2( Z ) Puc{(Z),a \M) = 
1. Hereafter, for notational simplification, the sum 



(Ei=i Z i = M (Z)). The standard deviation is calcu- sign of Eq. (|) will be written as J2 



lated according to the measure 



{N\M,(Z),<j} 



and 



\ 



1 M 

m E(^» 5 



(2) 



The equations obtained with the un-biased estima- 
tor of the standard deviation, used in Ref. Jl^] are 
listed in Appendix. The choice of the expression of 
the standard deviation will not have any influence on 
the shape of the correlation functions nor on the con- 
clusions of this study. It can be shown that the prob- 



the other sum signs will be formed according to the 
same logic. 

The extension of this formula of the denominator 
to samples containing a variable number of IMF is 
useful when the experimental statistics is reduced. 
It is expressed straightforwardly as P uc ((Z),a) = 
J2 M p M{M)P uc ((Z),a \M) (where P M is the mul- 
tiplicity probability distribution of the imf), i.e. : 



P uc ((Z),a) = J2 ^m(E^)(E^) ! 

{N|<Z),cr} Z Z 

Pz(Z\Ez' Nz-) Nz 



n 



N z \ 



(G) 



3 In the following, the variables in capitals will be rela- 
tive to the IMF, the IMF partitions will be noted N, the 
complete partitions n and the total multiplicity m. 
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B. Statistical error bars 

Let us recall that the correlation function is defined 
by: 



l + R((Z),a \M) = 



P((Z),<? \M) 
P uc ((Z),a \M) 



null hypothesis, that the peak has a height lower than 
that observed. Therefore, the higher the peak, the 
higher the significance. An under- estimation of the 
significance S can be obtained straightforwardly using 
the Schwarz inequality : 



(7) 



where the probability in the numerator is the num- 
ber of sample events including M IMF, with mean 
charge (Z) and standard deviation a, divided by the 
number of events with IMF multiplicity M . To first 
order, the sampling variance of a proportion applied 
to Eq. (|t]) gives the following error : 



A(l + R((Z),a \M)) 



y/P((Z),a \M) 



^N{M) P uc ((Z),a \M) 



(8) 



where N(M) is the number of events with M IMF in 
the data sample. The use of formula (||) reduces con- 
siderably the statistical error. In the case of a Monte- 
Carlo selection process, it would be, to the same or- 
der : 



A(l + R((Z),a |M)) = 



y/P((Z),a \M) 



y/N{M) P uc ((Z},a \M) 
P((Z),a \M) 



^/N UC (M) P uc ((Z),a |M)3/a ' 

where N UC (M) is the number of pseudo-events gen- 
erated by random selection for the calculation of the 
denominator. The last term can be very important 
in the presence of correlation peaks. The calculation 
of the error is crucial when the standard deviation is 
zero, on the one hand because it is these events that 
we are interested in, and, on the other hand, because 
the number of events of this type is often very small. 

In practical cases, the denominator may be evalu- 
ated with a very low uncertainty thanks to Eq. (|^) 
since only the - very low - statistical fluctuations on 
the charge spectrum alter the result. On the other 
hand the precision of the numerator depends strongly 
on the number of events in the considered sample. 
Furthermore, the error on the error bar (Eq. (||)) 
depends also on the number of events. Therefore it 
can be inaccurate. It would therefore be interesting 
to obtain an evaluation of the error bar using only 
the value of the denominator. This is possible using 
the so-called null hypothesis, i.e. that the correlation 
function is equal to unity (absence of correlation) . In 
the frame of this hypothesis the error bar is : 



A 1+R=1 (l + R({Z),a \M)) = 
1 

y/N(M) P UC ((Z),(T \M) 



(10) 



The significance of a positive correlation (of a peak) 
is defined as being the probability, in the frame of the 



S < 



N(M) P uc ((Z),<t | A/) 3 
(P((Z),a \M)-P nc ((Z),a \M)f 



(11) 



Exact calculations of the significance as well as ap- 
plications to experimental data will be presented in a 
forthcoming publication jlj]] . 



C. Case where all IMF have the same charge 

1. Numerator 

Since the spinodal decomposition peak is expected 
when all IMF have the same charge, we now con- 
sider the case where er = 0. For a fixed IMF mean 
charge, there is now only one IMF partition : Vi, Z\ = 
(Z). Thus, differences between the complete parti- 
tions with same (Z) are only due to the light frag- 
ments whose total charge is zn g ht = ^tot — M(Z). 



(9) 



2. Denominator 

When cr = 0, the probabilities given by Eq. (||) 
become : 

= P Z ((Z) \M) M 
if (Z) is integer and 
P uc ((Z),0 \M) I (Z) E [Z min ,z tot - (M - I) 

= o 

otherwise , 

(12) 

where Pz{Z\M) is the charge distribution for a 
given IMF multiplicity (£z=z min P{Z\M) = I). The 
mean charge being equal to the charge of each IMF 
implies that (Z) is always an integer. The probability 
that the standard deviation is zero is : 



Puc(<T = 0) = 



t -{M-l)Z„ 

E 



P Z {{Z)\M) M . (13) 



When the charge distribution of light IMF follows 
a power law or an exponential law (we will see that 
this is the case for the minimum information model), 
the denominator assumes very simple forms (respec- 
tively) : 



p uc ((z),o \m) oc (zy™, 

P uc ((Z),0 |M) oce- M < z > 



(14) 
(15) 
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3. Correlation function 

As indicated previously, the evaluation of the cor- 
relation function in the case of equal size IMF, is of 
considerable physical interest. Unfortunately it often 
corresponds to the ratio of two very small probabil- 
ities. If the number of events in the sample is too 
low, it is possible that no event corresponds to the 
given mean charge (the correlation function cannot 
be calculated) or that a very small number of events 
correspond (which can lead to a spurious peak). It is 
therefore important to determine, a priori, the mini- 
mum number of events necessary to obtain a reliable 
evaluation of the correlation function for a null stan- 
dard deviation. An evaluation of this number can be 
obtained making, once more, the hypothesis that the 
correlation function is unity. The probability to ob- 
tain an event in which the IMF charges are all equal 
to (Z) is P Z {(Z)\M) M (this quantity can be obtained 
precisely even with a reduced event sample) . The min- 
imum size N(M) of the sample must be therefore of 
one order of magnitude greater than P%{{Z)\M)~ M . 



III. DENOMINATOR CONDITIONED BY 
CHARGE CONSERVATION. 

The formation of the denominator as proposed by 
Moretto and collaborators (that we will continue to 
call the pseudo-event method though the result is ex- 
pressed by the algebraic formula (0)) has many ad- 
vantages : it is rigorous, it is simple to evaluate, it 
takes into account the efficiency of the detector, it 
uses only experimentally measured quantities and the 
resulting correlation function shows all charge correla- 
tions, whatever their origin. This latter advantage can 
become an inconvenience when one wishes to study 
correlations induced by only one physical cause. In 
the majority of cases, the main structure in the cor- 
relation function is due to the total charge conserva- 
tion law. We will see that this law introduces a large 
structur e, g reater than unity, close to (Z) = z tQt /M. 
In Ref. this structure was considered as a base- 
line on which was superimposed a peak due to the 
spinodal decomposition mechanism. 

In this section we will discuss two different propo- 
sitions for calculating the denominator taking into 
account the charge conservation (in order to remove 
the corresponding structure from the correlation func- 
tion). The first one consists in using partitions pro- 
vided by the minimum information model (all the par- 
titions of a given total charge have the same probabil- 
ity). We will show that the denominator constructed 
in this way presents a spurious peak at a = that 
can conceal a pos sible p hysical peak in the correlation 
function (Section [II A ). The second proposition con- 
sists in modifying the expression (H) of the algebraic 



calculation of the denominator in order to introduce, 
in an exact way, the influence of charge conservation 
with the consequence that charge conservation influ- 
ences both numerator end denominator. 



A. Minimum information model 

1. Introduction 

In this model, all partitions have the same proba- 
bility : P(n) = l/iV(z to t) where N{zt Q t) is the total 
number of partitions for total charge z to t- This re- 
sult is obtained by application of the minimum infor- 
mation principle (or maximal entropy), information 
being defined as : 



I = J2 p ( n ) lnP(n). 



(16) 



Setting the derivative of / equal to zero, under the 
single constraint of charge conservation, one obtains 
that all probabilities are equal. The total number of 
charge partitions of a charge Ztot nucleus is approx- 
imately given by the Ramanujan-Hardy formula |is}| 
whose leading term is : 



N(z tot ) 



exp 



(V 



3 



4z t otV3 



(17) 



The number of partitions increases therefore very 
rapidly with the charge. Thus, studies of large sys- 
tems, by systematic generation of all partitions, are 
not possible. The calculation of the number of IMF 
partitions (all fragments have a charge greater than a 
certain limit) and of light fragments (fragments with 
charge lower than a certain limit) is exposed in the 
companion article p9| . Some examples of applications 
of the minimum information model (possibly modified 
by combinatorial factors) to nuclear fragmentation are 
given in Refs. MM 



2. Case where all IMF have the same charge 

In this subsection, the correlation function for the 
minimum information model in the case a = will 
be calculated exactly. It will be shown that this func- 
tion presents a combinatorial peak due to an intrinsic 
feature of the model, namely the non-ordering of the 
charges. We will introduce an alternative model in 
which this effect is corrected. 

a. Numerator. The numerator P((Z),0 \M) of 
the correlation function is calculated as the number 
of partitions with M IMF, each of charge (Z), divided 
by the total number of partitions with M IMF. The 
charges of the IMF being fixed, the number of parti- 
tions will be equal to the number of ways to divide the 



5 



remaining charge zn g ht = z to t — M{Z) into light frag- 
ments (i.e. fragments with charge less than or equal 
to z max = Z m i n — 1). This constrained number of 
partitions will be noted 2max -/V(zii g ht). Similarly, the 
number of partitions of charge z into M fragments 
with charge greater or equal to Z m i n will be noted 
z min N(z. M). These numbers can be calculated ex- 
actly jlSfl. With our notation, the numerator reads : 



P((Z),0 \M) 



hi. 



{zy 



N{M{Z)',M) *™*N(z to t ~ M(Z)>) 



(18) 



b. Denominator. We have seen (Eq. (P-2|)) that 
the denominator is written as Pz((Z) \M) M when the 
standard deviation is zero. The conditional probabil- 
ity of (Z) given M is the number of partitions with 
M IMF weighted by the proportion of charges (Z) that 
they contain, divided by the number of partitions con- 
taining M IMF, so that : 

P uc ((Z),0 |M) 
J2 Z (z)Jz miD J2{n\(z),m} ~ J m~ Zta °-*N(z tot - M(Z)) 
E(zy=z min z min N(M(zy,M) ^N(z tot - M(Z)<) 

in which the sum over all partitions containing M 
IMF has been written : 



ztot/M 

E E = £■ 

(Z)=Z Iaia {N\{Z),M} Zi 



R k /(M-k+l) 

E -E> ( 2 °) 



with Zq = Z min and Rk = z tot - J2i=i The 
charges of the IMF are noted Zi and are written in 
increasing order. 

c. Correlation function. The charge correlation 
function is thus given by : 



1 + R((Z),0 \M) 



N(z tot - M(Z)) 



z mln N{M{Z)',M) z —N(z tot - M{Z)') 

(zy 

Ztat/M 



Af-1 



M 




io- 



A 
N 

V 



10 



M 




A 
N 

y_ 

+ 10 



(19) 
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FIG. 1. Numerator (upper row), denominator (median 
row) and correlation function (lower row) for total charges 
56 (left column) and 104 (right column). The bold lines 
correspond to IMF with minimum charge 5, the thin lines to 
IMF with minimum charge 3. The dotted lines correspond 
to the events with 5 IMF and the full lines to those with 3 

IMF. 

We observe that : 

• The behavior of the correlation function de- 
pends only weakly on the size of the system. 
This result has been observed experimentally : 
similar correlation functions have been observed 
for very different systems in central collisions 
H2l. 



N. 



(z) 



<N(z tot -M(Z))\ (21) 



(Z)=Z min {N|(Z),M} 



The latter result is, of course, free of error since it 
results from the numbering of all possible partitions. 
Results for two total charges, two multiplicities and 
two definitions of the IMF are presented in Fig. |l|. 



The denominators are exponentially decreasing. 
This is due to the fact that the charge distribu- 
tions are also exponentially decreasing between 
and ztot - (M - l)Z min (Eq. «). 

The abscissa of the maximum of the correlation 
function is a few units lower than z tot /M, this 
property can be used to provide an experimental 
determination of the total charge of the compos- 
ite nucleus, after pre-equilibrium emission, mea- 
suring only the IMF. 

The amplitude of the correlation function in- 
creases strongly as the minimum charge of the 
IMF diminishes and their multiplicity increases. 



Fig. g presents the correlation function obtained for 
different multiplicity 3 systems (left column) and for 
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the total charge 79 and different multiplicities (right 
column) . 



A 
N 
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<z> 



FIG. 2. Charge correlation function, for a null stan- 
dard deviation, obtained by formula ( ^l| ) in the framework 
of the minimum information model. Heavy traces corre- 
spond to IMF with charges greater than 5 and light traces 
to IMF with charges greater than 3. The correlation func- 
tions for an IMF multiplicity of 3 are presented on the left 
hand column in which the total charges of systems are, 
from top to bottom, 79, 104 and 400. The three graphs on 
the right hand column correspond to the total charge 79 
and the IMF multiplicity runs from 4 (top) to 6 (bottom). 
The horizontal lines are placed at 1 + R = Ml. 

The left hand column shows that the correlation 
functions are nomothetic for a given multiplicity. 
On each figure a horizontal line has been placed at 
1 + R({Z),0\M) = Ml. One oberves that the maxi- 
mum of the correlation function is always of the order 
of Ml. This property is due to the multinomial factors 
of the denominator (Eq. (|^)). When the standard de- 
viation is null, the product \\ z Nzl is identically equal 
to Ml, it factorises therefore in the numerator of the 
correlation function. When the standard deviation in- 
creases, this product decreases rapidly down to 1 when 
all charges are different. To illustrate this point, let us 
consider two very similar partitions of the total charge 
21 into 3 IMF : {7,7,7} and {6,7,8}. The numerator 
of the correlation function for both partitions is the 
same : l/iV(21,3) (= 1/249). The denominator is, 
in the former case, P(7|3) 3 (= 6.6 10~ 4 ) and, in the 
latter case, 3! P(6|3)P(7|3)P(8|3) (= 3! 6.04 1CT 4 ). 



Therefore the probability product is almost the same 
in both cases. Besides, the charge conservation con- 
straint is weak when the IMF have small charges hence 
1 + P({6, 7,8}) w 1 (one finds 1.1). Consequently 
1 + R({7, 7, 7}) w 3! (one finds 6.1). 

It can be conjectured that this effect would occur 
for all non-ordered fragment models. If one multi- 
plies every partition by a factor equal to the number 
of charge permutations (Ml/Y\ z Nzl), the peak dis- 
appears. On the other hand probabilities for other 
sigma values arc little modified because the Nz are 
then almost always equal to or 1 (i.e. Ylz ^z* = !)• 
This latter model will be referred to as the ordered 
minimum information model. 



B. Algebraic calculation of the denominator with 
charge conservation 

1. Method 

The goal of this subsection is to introduce the ex- 
act method for the evaluation of the denominator that 
eliminates the effects due to charge conservation from 
the correlation function. This denominator is ob- 
tained by an extension of the formula (H) to the whole 
charge (i.e. including light fragments). In a first step, 
we suppose that there is no correlation at all between 
charges. This means that each charge is described 
by a probability lntr P z (z) (that will be referred to as 
intrinsic probability of the charge). The conditional 
probability of a partition n (including the IMF and the 
light particles) with total multiplicity m is given then 
by the multinomial formula : 



P(n|m) = ml TT ■ 



(22) 



These conditional probabilities obey the normaliza- 
tion condition X)n^( n l m ) = ^ one introduces 
the constraint of total charge conservation, partition 
constrained conditional probabilities are given by (an 
index cc will be applied to probabilities constrained 
solely by the charge conservation) : 



P cc (n\m) = k(m) ml z -\-^- — S 



with 



k 1 {m) — ml i J • 



mtr P z (z)" 



(23) 



(24) 



{n|m,2 tot } 2 



On the other hand, the multiplicity probability dis- 
tribution is given by : 



7 



P m (m) 



e n 

{n\m,z tot } z 



l 'Pz(z) r 



with 



«- x = e (E^n- 

{n|z tot } z z 



^zW 



(25) 



(26) 



Finally, the partition probabilities are given by 



intr p / v \n z 



(27) 

These probabilities contain all the information rel- 
ative to the charges and to their correlations. For ex- 
ample, (observed) charge and conditional charge prob- 
ability distributions are given respectively by : 



E n^ P cc {n) 

Pz{za} = ^„(i>*)^w 

= E{n|z tot } n ^o (Ez n z) ! II; 

Pz(zo\m) = p-^Y,{n\m} ^ P cc(n) 



intI P,(z)"" 



(28) 



Because of correlations, the intrinsic probabilities 
arc not equal to the probabilities to observe the 
charges (the conservation constraint favours small 
charges, see Fig. ||b). Equality between intrinsic prob- 
abilities and observed probabilities is valid only for an 
infinite system. 

The intrinsic probabilities are quantities which are 
not directly measurable, they must be calculated by 
inversion of Eq. (|2^) where the P cc (n) are the mea- 
sured frequencies of the partitions. The set of Eqs. 
( ^8| ) constitutes an under-determined system. It is 
thus not possible to obtain an unique solution. 

However, inversion of (^7|) is possible if the non triv- 
ial correlations between charges in the studied sample 
are weak. When the intrinsic probabilities are deter- 
mined, the probabilities of the denominator can be 
calculated by summing the complete partition prob- 
abilities having the same IMF mean charge and stan- 
dard deviation : 



P cc ((Z),a\M) 



Pm(M) 



E P -(n). (29) 



{n\M,{Z),a} 



Finally we can write that the denominator is con- 
structed using the intrinsic probabilities : 



E 



{n\z tot ,M,{Z),<y} 



(E, 



P cc ((Z),a\M) 

-z)inz^tF 



E< 



t,M} 



(Ez«z)!n 2 



tr -P,(z)"z 



(30) 



This new denominator takes explicitly into ac- 
count the charge conservation. Structures observed 
in the corresponding correlation function will neces- 
sarily arose from other causes. 



2. Monte-Carlo generation of events without non trivial 
charge correlations 

A sample of events without non trivial correla- 
tions may be synthesized by the following procedure. 
Charges Zi are selected randomly according to their 
intrinsic probabilities until Zi > z to t ■ The event is 
preserved only if J2i z i — z tot ■ The resulting charge 
spectra are those given by Eqs. (p8|). The alternative 
procedure that would consist in randomly selecting 
M — 1 charges, and deducing the last charge using 
charge conservation would introduce a bias since the 
distribution of the last charge would be different from 
the preceding ones. 



3. Combinatorial independent emission models 

We stated previously that models can be charac- 
terized by the term "independent emission" if they 
contain no correlations, other than those induced by 
the conservation of the total charge. This definition 
implies that partition probabilities generated by such 
models can be written in the form of Eq. (p7|). We 
give here two examples of independent emission mod- 
els. 

a. Ordered minimum information model. The 
model produced by weighting partitions by a factor 
m VIlz riz ' ^ s an independent emission model. One 
notes that this weighting is the same as the one of 
Eq. @ if the product ]J Z intI P z (z) n * is constant. 
This condition is fulfilled if and only if lntT P z (z) = a z 
(the intrinsic probability product is a Ztot ). Thus, the 
normalization of the probabilities is Ez flZ = For 
a sufficiently large total charge (in practice, superior 
to 10), the normalization condition implies a = 1/2. 
The model in which every partition is weighted by the 
number of permutations of its charges (to!/ Ilz n J) 
is therefore an independent emission model with in- 
trinsic probabilities lntT P z (z) = 2~ z . The resulting 
charge distributions are almost exponentially decreas- 
ing (Fig. Conversely, the minimum information 
model is not an independent emission model (it is not 
possible to find a set of intrinsic probabilities such that 
all partition probabilities would be the same). 

b. Model of charge equi-probability. The mini- 
mum information model is taken as implying the equi- 
probability of the partitions. In this context it is also 
interesting to see results given by another elementary 
model obtained by assuming intrinsic equi-probability 
of the charges. The resulting charge distribution is of 
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course not uniform because of the constraint of con- 
servation of the total charge. 

The Eq. ©, for intI P z (z) = l/z tot , gives : 



F cc (n) = a(£,n z )!n z H# 



^ rn 



(31) 



In other words, if one generates partitions by im- 
posing charge conservation and if all charges have the 
same intrinsic probability to be selectioned, then the 
partition weight is m!/(z[£ t fl z n z \). Examples of ob- 
served charge spectra are given for different values 
of the total charge in Fig. ||b). These spectra are 
very different from that of the intrinsic probabilities. 
Indeed one observes an almost exponential decrease. 
Only the last charge has a probability which is not 
consistent with this tendency. This behaviour is due 
to the fact that the charge Ztot — 1 can appear ac- 
companied only by a charge 1, whereas the charge 
z to t is never rejected. More precisely, Eq. ( |28| ) gives 
Pz{ztot) = (^tot/2) P z (z to t - 1). One notices that 
the greater the total charge the less the influence of 
the constraint of conservation and consequently the 
smaller the slope of the exponential. 
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FIG. 3. Charge distributions obtained in the or- 
dered minimum information model (a) and in the charge 
equi-probability model (b) constrained by total charge 
conservation. Circles correspond to a total charge of 20, 
squares to z to t = 56 and points to Ztot = 104. The full 
lines are the intrinsic probability distributions. 



4- Multiplicity constrained independent emission model 

We now consider the model for which probabilities 
of partitions with fixed multiplicity are given by Eq. 
fl22|), but the multiplicity probability distribution is 
not given by the combinatorial Eq. (|2^). The mul- 
tiplicity probability distribution P m (m) is imposed a 
priori. This model has been studied extensively by 
A.J. Cole and collaborators (24). The main difference 
between the quoted articles and the study presented in 
this paragraph resides in the interpretation of quan- 
tities noted X z by Cole et al. and, here, lntr P 2 as, 
respectively, adjustable parameters and intrinsic prob- 
abilities. However, the X z parameters being defined 
to a factor a z , it is always possible to normalize them. 
Indeed, the product T\ z {a z X z ) nz is equal to the prod- 
uct Y[ z X?* to within a constant (a Ztot ). Eq. @ and 
the normalizations give : 



P cc (n) 



P m (m) 



n 



Zj{n'|m,z t ot} -112 

iatr P z (z) n * 



with 



m = > n 



(32) 



(33) 



In the same way as for charges, one can introduce an 
intrinsic probability for multiplicities ( mtv P m )- The 
partition probability is then written as : 

intr p / 7 \n z 

P cc (n) = k l ^P m (m) m \ J] J ) * ~ znz 

2 

(34) 

with : 



intr p / -,\n z 



{n|m,2:tot} z 



(35) 



One deduces the relation between observed proba- 
bilities of multiplicities and those of the intrinsic prob- 
abilities : 



P m (m) = mtr P m (m) 
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{n|m,2 t ot} 



il 



a 



»"P ; (z)"i 



Intrinsic probabilities and observed probabilities 
possess the same distributions in the limit of an in- 
finite system size. 



IV. APPLICATIONS 



prediction which is verified by experimental observa- 
(36)tion. The first two factors can be considered therefore 
to be constant for a given multiplicity. The volume 
of the nucleus at the time of fragmentation, Vp, is 
supposed to be independent of the partition (or, ac- 
cording to the version of the model, dependent only 
on the multiplicity). If we disregard the conservation 
of the mass, the emission probability of a fragment 



A. Introduction 

In this section we show the charge correlation ob- 
tained, for several nuclear decay models, using the 
denominator given by the independent emission hy- 
pothesis. The first step of the procedure consists in 
determining the intrinsic probabilities of the charges 
for each model sample. These probabilities are ob- 
tained by a recursive procedure of minimization of 
the x 2 between probabilities of partitions in the syn- 
thesized sample and those given by Eq. (|2^). The 
convergence of the procedure is possible only if non 
trivial correlations between charges are weak. The 
minimum x 2 is therefore an indication of the strength 
of these correlations. The second step, calculation of 
the denominator by the method presented in Section 
III B 1| pOl) , would also not be possible in the presence 



of strong correlations. We will see that this condition 
is fulfilled by the three models that we are going to 
study. Results of the application of this procedure to 
the experimental events will be presented in forthcom- 
ing articles Jl7|,p5[. 



B. The Copenhagen model 



is proportional to a 3 / 2 : P z {z) = aa?/ 2 : so that the 
product Y\ z Pz{z) 7lz = a m Y[ z cl 3t1z ^ 2 involves a new 
factor depending only on the multiplicity. The equa- 
tion can therefore be rc-writtcn as : 



1L nJ. Zz ">»'L ! 



(38) 



so that we recover the same expression as that ob- 
tained for the independent emission model. We thus 
expect that the correlation function is everywhere 
equal to 1 and if the denominator is calculated from 
pseudo-events, we expect that the shape of the corre- 
lation function is determined by charge conservation. 
However, this conclusion is based on a simplified form 
of the model. The correlation function may thus ex- 
hibit weak modulations. Moreover, one does not ex- 
pect a peak for small values of the standard deviation. 

Using the SMM code, 35 million events have been 
generated for the 138 Ba nucleus excited to 5 MeV 
per nucleon. We built, from this sample, the charge 
correlation function for the IMF. The Copenhagen 
model produces results (almost) consistent with in- 
dependent emission as shown in Fig. ^. Discrepancies 
can be explained notably by the fact that the hot frag- 
ments produced during the multifragmentation phase 
described by Eq. (|37j), thereafter decay by light par- 
ticle emission. 



The Copenhagen model |]]2(| is a hot liquid drop 
model that describes the multifragmentation of the 
nucleus as an instantaneous statistical mechanism. 
The probability of a partition in mass and in charge 
n a ,z for an excitation energy E* is given by : 



P(n\E*) = P!(E*) 



Vf 



x<5„ 



n 

,a,z 

4, E» 



7 3n a>i /2 



! 



(37) 



The two Kronecker symbols account for the conser- 
vation of mass and of charge . The de Br oglie wave- 
length of the nucleon At = \/2ttH/ (m n T) depends on 
the temperature which is roughly constant at fixed 
excitation energy. The density of states correspond- 
ing to the internal excitation energy of the fragments 
Pj(E*) is also constant at a given excitation energy. 
In this model, as in comparable models, the multiplic- 
ity is correlated linearly with the excitation energy, a 




5 10 15 20 25 30 



FIG. 4. Comparison of charge spectra for several mul- 
tiplicities for SMM (circles) and for the independent emis- 
sion hypothesis (lines) . The dashed line corresponds to the 
intrinsic probabilities. The multiplicity probability distri- 
butions, P m (m), given by the two models are also shown 
in the figure. 
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The charge correlation functions calculated with the 
denominators of formulas (||) and (|3^) are presented in 
the Figs. |5| and |[ In the first case, the main structures 
are due to the conservation of the charge, in the sec- 
ond, the correlation function is practically flat. Non 
trivial correlations between charges are therefore very 
weak for this model and no preferential fragmentation 
into equal charge is observed. The - small - modula- 
tions of the correlation function are due to physical 
causes and to statistical fluctuations. 




FIG. 5. Upper and lower figures present, in order, the 
numerator, denominator and [{Z),a] correlation function 
for the Copenhagen model. The denominator is calculated 
by the method of pseudo-events (formula (^)). 




FIG. 6. Denominator (upper figure) and correlation 
function (lower figure) of [(Z),cr] for the Copenhagen 
model (the numerator is presented at the top of the pre- 
vious figure). The denominator is calculated using the 
independent emission hypothesis, via formula (po[). 

C. Percolation 

The same study has been carried out on a sample 
obtained with a percolation code jL6|. A sample of 
10 8 events has been generated using a 3D percolation 
program on a simple cubic 4x4x4 periodic frame, for 
a bond breaking probability of 70 %. The result of 
the \ 2 minimization process is presented Fig. |?] in 
which are compared charge distributions for various 
multiplicities in percolation and those given by the 
formula (^3|) (the intrinsic probabilities are indicated 
by a dotted line). The very good agreement indicates 
that the correlations between charges are very weak 
in this model. 
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FIG. 7. Percolation (same conventions as in Fig. y.) 
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Figs. H and || present percolation correlation func- 
tions obtained using denominators given respectively 
by the pseudo-event method (Eq. (||)) and by the 
independent emission hypothesis (Eq. (j3p|)). These 
figures were constructed for all possible values of the 
number of IMF. In the first case, structures are almost 
entirely due to the conservation of the charge (to each 
IMF multiplicity corresponds an edge line). When the 
denominator is calculated using the intrinsic proba- 
bilities, the correlation function is flat and equal to 1 
(Fig. |^). The small peaks on the sides of the correla- 
tion function are due to the statistical fluctuations. 

The model of percolation can therefore be assimi- 
lated to knowledge of the intrinsic probabilities. 
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FIG. 8. Percolation (same conventions as in Fig. 




FIG. 9. Percolation (same conventions as in Fig. ^.) 

Fig. [h] presents the correlation function of the per- 
colation calculation when the denominator is calcu- 
lated with the minimum information model (upper 
figure) and with the ordered minimum information 
model (lower figure) for same total charge (ztot = 64). 
In both cases, the correlation function presents large 
structures which are not easily interpreted. 




FIG. 10. Correlation function of [(Z),cr] for the per- 
colation process (the numerator is presented at the top of 
Fig. ^). Denominators are given by the models of the 
minimum information (upper figure) and of the ordered 
minimum information (lower figure). 
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D. The Brownian One Body dynamical model 

The last model to be studied in this work is char- 
acterized by non trivial correlations, even though, as 
motioned in the introduction, they are partly masked 
by different processes. The preferential decomposition 
of the system into almost equal charges (in the range 
[10,20]), which characterizes the BOB model, is not 
visible, for example, in the inclusive charge spectra 
(Fig. 0). 

Our simulated sample was obtained via a four step 
process |27j . The collision entrance channel has been 
simulated using a one-body semi-classical microscopic 
calculation of the bnv type |^8| for the system 129 Xe 
+ 119 Sn at 32 A MeV. This calculation shows that, 
for the most central collisions, a compressed single 
source is formed, after a weak pre-equilibrium emis- 
sion, within 40 fm/c. The decompression phase, the 
entrance into the spinodal zone and the formation 
of the fragments are then followed by the BOB code 
which simulates Boltzmann-Langevin density fluctu- 
ations [^9| and the evolution of the system density 
submitted to spinodal instabilities. In a third step, 
fragments are formed using an algorithm which re- 
groups contiguous cells which numbers of test parti- 
cles (40 test particles are used to simulate a nucleon) 
is greater than a given threshold. The resulting nu- 
clei are hot, their statistical decay and their Coulomb 
expansion are, in a last step, simulated by the SIMON 
code The event samples are eventually filtered 

using the indra response function po| . 

The regrouping of pseudo-particles generated by the 
model is possible only for fragments of charge superior 
or equal to 5. The light fragments are not known. This 
difficulty must be taken into account in the routine of 
intrinsic probability optimization : the probability of 
a partition of M IMF is the sum of probabilities of all 
partitions containing these M IMF (solely) together 
with the corresponding light particles : 



P(N)=a Yl (£« ! I1 

{n|2 t ot,N} Z Z 



intr 



(39) 



The calculation of these probabilities can be acceler- 
ated considerably by noticing that they can be written 
in the form : 

x ( £ c£f +m k(m, z tot - £ ZN Z )\ 

\ m Z / 



(40) 



The last factor depends only on the sum of charges 
of the IMF and on their multiplicity (k is given by Eq. 
(pi|)). The result of the fit of the intrinsic probabilities 
is given in Fig. O. 
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FIG. 11. Model Brownian One Body (same conven- 
tions as in Fig. |i|.) 

One notices that, in spite of the absence of light 
particle in the sample, the intrinsic probabilities have 
a realistic distribution for all charges. The resulting 
correlation function is presented in Fig. 



12. The par- 



tition probabilities are correctly reproduced by the 
independent emission hypothesis, hence the correla- 
tion function is practically flat. However, in contrast 
with the previous models, it includes strong correla- 
tion peaks near a = 0, and, to a lesser extent, for the 
maximal values of a for given (Z). These latter peaks 
(as well as those corresponding to a = 0, (Z) < 9) 
have a low significance, so they can only be due to the 
statistical fluctuations. Peaks at a = 0, (Z) > 10, on 
the other hand, arc meaningful. They signify the spin- 
odal decomposition of nucleus produced by the bob 
code. 
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FIG. 12. Brownian One Body model correlation func- 
tion (same conventions as in Fig. ^|) 



V. CONCLUSIONS 

The main goal of this article was to present a new 
method for the evaluation of the non trivial correla- 
tions between the fragment sizes of a finite size system. 
The conclusions are the following : 

• The Monte-Carlo calculation of the denominator 
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proposed by Moretto and collaborators can be 
replaced by a fast algebraic calculation which is 
equivalent to the selection of an infinite number 
of pseudo-events (||). This calculation results in 
decreased error bars in the correlation function 
notably those associated with correlation peaks 
(||) . The extension of the calculation of the cor- 
relation function to samples including a variable 
number of IMF presents no particular difficulty 




Correlation functions thus obtained for the dif- 
ferent models studied in this article (SMM, per- 
colation, minimum information, bob) all possess 
one maximum for (Z) smaller by a few units 
than ztot/M. This property leads to the evalua- 
tion of the size of the composite nuclei for which 
decays have been observed experimentally. 

This type of denominator possesses many advan- 
tages. However, correlations induced by charge 
conservation are always important. They may 
conceal other less trivial correlations, or distort 
the evaluation of their amplitude. It would, 
therefore, be useful to define a conventional 
method for calculation of the denominator in- 
cluding effects induced by the charge conserva- 
tion. 

It has been proposed to evaluate the denom- 
inator using the minimum information model 
(all possible partitions of a given total charge 
has the same probability). This model incorpo- 
rates charge conservation but possesses a purely 
combinatorial correlation peak at a = so 
that there is a risk of concealing a physical 
peak present in the data sample (this effect 
can be corrected by the weighting of the par- 
tition probability by the number of its permuta- 
tions : Tn\/ Y\ z n z ). Furthermore, the numera- 
tor (physical sample) and denominator (given by 
this model) may correspond to distinct charge 
and multiplicity distributions. Finally, the cor- 
relation function presents numerous structures 
which are difficult to interpret. 

The two previous conclusions lead us to propose 
a new calculation of the denominator, the goal 
being to replicate all features of partitions of 
the numerator excluding intra-event correlations 
due to other reasons than charge conservation. 
In the case where these non trivial correlations 
are weak, this goal is reached exactly using the 
independent emission hypothesis constrained by 
the conservation of the charge. Probabilities of 
partitions are given by the formula ( p7| ) which 
is based on the specification of intrinsic proba- 
bilities for each charge. These values represent 
probabilities for a charge to be observed if the 



constraint of charge conservation played no role. 
The intrinsic probabilities are not observables, 
so that they must be searched for by a proce- 
dure of minimization between probabilities of 
partitions in the data sample and those given 
by the formula (p7|). If the resulting x 2 1S low, 
it means that the studied sample is essentially 
composed of events corresponding to indepen- 
dent emission. The partition correlation func- 
tion (i.e. the set of ratios of the probabilities of 
the sample divided by the probabilities given by 
Eq. (p7f)) must then be always near unity ex- 
cept possibly for a reduced number of partitions 
corresponding to the non trivial correlations. In 
this work, only the [(Z), a] correlation has been 
studied, but the same procedure can apply to 
any type of correlation. 

• The proposed method has been applied to three 
models of nuclear multifragmentation. It has 
been shown that all three models correspond 
to almost independent emission. The first two 
(percolation and the SMM statistical multifrag- 
mentation code) result in correlation functions 
everywhere equal to 1 (to within 10 %). The 
code bob, on the other hand, exhibits a flat cor- 
relation function everywhere except at a = 0. 
These correlation peaks are due to the mech- 
anism of spinodal decomposition that favours 
partitions which include IMF of the same charge. 
These results legitimate the use of the charge 
correlation functions method for the experimen- 
tal search for spinodal decomposition. 

In forthcoming articles we will study problems re- 
lated to the application of this method to experimen- 
tal event samples (superposition of sources , distribu- 
tion of total charge, pre-equilibrium emission, experi- 
mental efficiency, calculation of the significance of the 
result) and we will present results obtained by the in- 
dra collaboration for heavy ion central collisions near 
the Fermi energy. 

I wish to thank I.N. Mishustin who suggested this 
work, M.F. Rivet, B. Borderie and M. Parlog for nu- 
merous and fruitful discussions and A.J. Cole for pre- 
cious advice concerning the manuscript. 



APPENDIX: EQUATIONS RESULTING FROM 
THE UTILIZATION OF THE NON-BIAISED 
ESTIMATOR OF THE STANDARD 
DEVIATION 

In this article, we used the usual definition of the 
standard deviation (||). The authors of Ref. (l3) pre- 
ferred to use the non-biased estimator (in this sense 
that the mean of its sampling function is equal to the 
real value). In writing this article, a being used as a 
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measure and not as an evaluation of the standard de- 
viation of an unknown distribution we restricted our- 
selves to its usual definition. In this Appendix we 
give the equations which result from the use of the 
non-biased value, er n b : 



(2) 



Cub 



Z i Zm - 1 

xP z (Z 1 \M)...P z (Z M - 1 \M) 
xPz({Z}'- 



(Al) 



M(M-2) ,„ , x 
Ma 2 nh - M _ 1 < \M) 



Mai 



M(M 



< 2 b \M) 



P uc ((Z),a nh |M) = £\..£ 

Z i Zm 

xP z (Z!\M) . . . Pz(Z M -i\M) P z {Z m \ M) 
xS. 
xS 



Z M ,M(Z) + (M-1)(Z)> 



z M< (zy±^M„i b -M^i <b 
d) -^Puo((Z),a nh \M) = M\ 

e n 

N Z 

2_^ z n z=m 

^ ? Z 2 N Z = M (Z) 2 + (M -l)^ b 



The other equations are not modified. 



P z {Z\M) Nz 
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(A3) 
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